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ABSTRACT 

The apparatus of correlation gamma function (r*(r)) is used to analyze 
volume-limited samples from the DR4 Main Galaxy Sample of the SDSS survey 
with the aim of determining the characteristic scales of galaxy clustering. Up to 
20/i~ 1 Mpc (Hq = 65 km s _1 Mpc _1 ), the distribution of galaxies is described by 
a power-law density- distance dependence, T*(r) oc r~ 7 , with an index 7 « 1.0. A 
change in the state of clustering (a significant deviation from the power law) was 
found on a scale of 20 — 25/i _1 Mpc. The distribution of SDSS galaxies becomes 
homogeneous (7 ~ 0) from a scale of ~ 60/i^ 1 Mpc. The dependence of 7 on the 
luminosity of galaxies in volume-limited samples was obtained. The power-law 
index 7 increases with decreasing absolute magnitude of sample galaxies M abs . 
At M abs ~ —21. 4, which corresponds to the characteristic value M* of the SDSS 
luminosity function, this dependence exhibits a break followed by a more rapid 
increase in 7. 

Subject headings: galaxies, groups and clusters of galaxies, intergalactic gas, 
large-scale structure. 
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1. Introduction 

The large body of observational evidence, in particular, the homogeneity and isotropy 
of the blackbody cosmic microwave background radiation, suggests that the distribution of 
visible matter is homogeneous on fairly large scales. The homogeneity of the matter distribu- 
tion in the Universe is postulated by the standard model of cosmological evolution (Peebles 
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1980). At the same time, the large-scale structuring of galaxies is traceable on scales larger 
than 300 Mpc (the Sloan Great Wall structure of of the SDSS survey). On small scales, 
the integrated parameters of the distribution of galaxies and systems of galaxies exhibit a 
power-law decrease in the density of objects as the volume under consideration increases. 
This behavior of the density is occasionally interpreted as fractality, self-similarity of struc- 
tures in a certain range of scales. Both the extent of fractal structures and the possibility 
of describing the observed distribution in terms of fractal sets (Davis 1997; Sylos Labini et 
al. 1998; McCauley 2002; Barishev and Teerikorpi 2005) are being debated. If fractality 
were found on significant scales, then its formation would have to be explained theoretically 
and reproduced in numerical simulations of cosmological evolution. A large extent of frac- 
tal structures (of the order of several hundred Mpc) would pose serious challenges to the 
standard theory of galaxy formation. 

To determine the characteristic scales of galaxy clustering, we use the conditional density 
function or the correlation gamma function (Coleman and Pietronero 1992). Based on a 
sample of rich Abell clusters, Tikhonov et al. (2000) determined the scale, ~ 100/i _1 Mpc 
(for the Hubble constant Hq = 100 km s _1 Mpc _1 ), from which the distribution clusters 
becomes homogeneous using this method. Based on galaxies from the CfA2 and SSRS2 
surveys and clusters from the APM survey, they established that the regions in which the 
density-distance dependence is well fitted by a single power law are limited by a scale of 
~ 30/i _1 Mpc. To test these results, in particular, to determine the homogeneity scale in 
the distribution of galaxies, we use data from the SDSS survey (York et al. 2000), which is 
currently the largest survey in sky coverage and depth. Based on the Luminous Red Galaxy 
Sample (Eisenstein et al. 2001) of SDSS galaxies, Hogg et al. (2005) found a homogeneity 
scale of ~ 70/i _1 Mpc (for Hq = 100 km s _1 Mpc _1 ). In this paper, we analyze a sample 
of galaxies from another part of the SDSS survey, the Main Galaxy Sample (Strauss et al. 
2002). 

Another principal objective of the correlation analysis is to describe the clustering of 
galaxies on small scales. The galaxy clustering parameters are known to depend on the lumi- 
nosity, morphological type, and color of galaxies (see Zehavi (2005) and references therein). 
Bright galaxies exhibit a stronger clustering than faint galaxies; the differences are more 
pronounced beginning from L*, the characteristic value of the Schechter luminosity function 
for galaxies. A detailed dependence of the clustering on galaxy properties is required to 
constrain the parameters of the theory of galaxy formation, in particular, to ascertain the 
initial conditions and to determine the formation conditions of these properties. This depen- 
dence was determined mainly using a two-point correlation function (see Zehavi (2005) and 
references therein). In this paper, we establish the dependence of the degree of clustering on 
galaxy luminosity using a correlation gamma function. 
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2. THE DATA 

When fully implemented, the SDSS (Sloan Digital Sky Survey) is expected to yield 
spectroscopic redshifts for ~ 10 6 galaxies and ~ 10 5 quasars based on photometric data 
from a sky region with an area of 10 4 sq. degrees in the Northern Galactic Hemisphere 
in five bands, u, g, r, i,and z, with a limiting magnitude of r = 22.5. The photometric 
data were used for a homogeneous selection of objects of various classes (galaxies, stars, and 
quasars) to be subsequently observed spectroscopically. Two types of galaxies were chosen 
to determine the redshifts from the list of objects classified as extended ones: galaxies with 
Petrosian magnitudes r < 17.77 and surface brightnesses higher than 24m/D" formed the 
Main Sample (the number of such objects in the final SDSS version is 900000) and the LRG 
(Luminous Red Galaxies) list includes galaxies with very red colors and r < 19.5 (the final 
SDSS version contains 100000 such objects). In this paper, we analyze data from the fourth 
SDSS data release DR4 (www.sdss.org) (849920 spectra, including 565715 galaxies and 76484 
quasars). 

When processing the DR4 data, we first selected a rectangular part from the region of 
spectroscopic sky coverage to make the allowance for the boundary conditions more conve- 
nient when calculating the gamma function and to ensure completeness of the sample. In 
the coordinate system of the survey (A, rj) (with the poles at 0:2000 — 95°,5 2 ooo — 0° and 
«2ooo = 95°,(5 20 oo = 0° 0-2000 = 275°,5 2 ooo = 0°; the point A = 0°,r] = 0° corresponds to 
02000 = 185°,<5 2 ooo = 32.5°, rj increases northward), the selected region is —47° < A < 8°, 
9° < i] < 36°. This region contains a certain number of small gaps in the spectroscopic 
coverage of the sky area under consideration. Filling the corresponding regions with the 
density equal to the mean density of the galaxy distribution in the selected region showed 
these gaps to have virtually no effect on the result. 

The Main Galaxy Sample is an apparent-magnitude-limited survey, which determined 
the method for constructing a volume-limited sample, determining the constraint on the 
absolute r-band magnitude of sample galaxies Mu m = ri im — 25 — 5lg(R max (l + z max )) — K(z)), 
where r^ m = 17.77 was taken as the limiting r-band magnitude and K(z) is the K-correction. 
We chose the sample parameters z max = 0.16 and z m i n = 0.05 based on the z-distribution of 
galaxies (histogram) so as to select the region of required completeness in redshift and to be 
able to fit fairly large spheres into the derived geometrical boundaries. 

We constructed a volume-limited sample containing 18916 galaxies with M abs < —21.77. 
The mean K correction for SDSS galaxies in the form K(z) = 2.3537,2 2 + 0.5735^ - 0.18437 
(Hikage et al. 2005) was used to estimate the absolute magnitudes of galaxies. We recalcu- 
lated the metric distances from the redshifts using the Hubble constant H = 65 km s _1 Mpc _1 
and the density parameters .0, vacuum = 0.7, f2 = 0.3 (see, e.g., Hogg 1999). 
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3. 



THE METHOD 



In the classical method of gamma function, the behavior of the density with distance is 
calculated in spheres (integral T*) and spherical layers (differential T). Any sample objects 
the spheres around which are located completely within the specified geometrical boundaries 
of the sample can be taken as the sphere centers. 

The differential (r(r))and integral (r*(r)) gamma functions are defined by the following 
formulas: 



where n{f) — i ^ fifa — f) is the number density ris the working radius, N is the number of 



the centers of the spheres involved in the counts, N is the number of objects in the sample, 
and A is the thickness of the spherical layer (it is assumed to be small). 

The gamma function gives the density of the objects within the spherical layer r < r < 
r-Q + dr around each object, i.e., the density at a given distance from the sample object for the 
differential gamma function and the density in spheres of radius r$ centered on the sample 
objects for the integral gamma function (the sphere centers are not included in the counts, 
i.e., the density of the "neighbors" is measured). In the case of practical implementation, 
the volume of the spherical layer is calculated as | • pi ■ (r 3 — Tq). 

The counts are averaged. If part of a particular sphere goes beyond the sample bound- 
aries as the working radius ro increases, then this sphere is excluded from the counts. Thus, 
as the working radius increases, only the spheres with the centers that come increasingly 
close to the center of the volume covered by the sample are involved in the counts. The 
counts are stopped when the number of remaining spheres N sp < 10. The scale on which 
the gamma function is used is limited by the radius R s of the largest sphere centered on the 
sample object that fits into the geometrical boundaries of the sample. 

The result is presented on a logarithmic scale as the dependencies of log{T) on log(ro) 
and of log(T*) on log(r ). The slope 7 of the fitting straight line that is constructed using the 
selected region of change in log(ro) determines the correlation dimension of the distribution 
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(codimension) D corr = 3 — 7 (7 > 0). If the distribution is fractal, then D corr = D, where D 
is the fractal dimension of the distribution. 

Generally, a higher slope (a higher 7 corresponding to a lower dimension) means a 
greater mean decrease in density inside the volume and, consequently, a larger clustering of 
objects. The horizontal portions of the plot indicate that the distribution of objects in the 
sample is homogeneous on the corresponding scales (D corr = 3). 

Clearly, the slope in the logarithmic dependence (a good linear fit) is insufficient to 
assert that there is fractality on the corresponding scales (McCauley 1997; Tikhonov 2002). 
Nevertheless, the method of gamma function is an appropriate apparatus for determining 
the characteristic scales of galaxy clustering irrespective of whether the fractal interpretation 
is valid or not. 



4. CHARACTERISTIC SCALES OF THE GALAXY DISTRIBUTION 

Figure 1 shows the pattern of SDSS galaxy density variations with distance. We can 
reliably distinguish the main features of the galaxy distribution on various scales. Up to 
20/i" 1 Mpc, the gamma function is well fitted by a power law with the index (the slope in the 
logarithmic axes) 7 = 0.98 ±0.02. On a scale of (20 — 25/i~ 1 Mpc, significant deviations from 
a linear dependence begin — the gamma function exhibits a break; the change in the slope 
is sharper for T(r), which reflects the change in the state of clustering more accurately than 
T*(r), since the latter is more inertial. This suggests that the corresponding changes in the 
galaxy distribution also occur on a scale of 20/zT 1 Mpc abruptly. This effect may be related 
to the sizes of superclusters. The transition regime of clustering follows next (the density 
decreases with distance, but more slowly than on the power-law portion) and the distribution 
of SDSS galaxies becomes homogeneous beginning from a scale of 50 — 60/i _1 Mpc. The depth 
of the sample allows spheres up to 130/zT 1 Mpc (the scale R s ) to be fitted in its volume and the 
extent of the "homogeneous" portion of the density variation with distance is considerable. 

It follows from our analysis of volume-limited samples with the angle boundaries —47° < 
A < 8° and 9° < i] < 36°, different constraints on the radial coordinate (z ma x), and, accord- 
ingly, different upper limits on the absolute magnitudes of the galaxies included in the sample 
(up to the sample with the largest sphere that can be fitted in its boundaries, ~ 100/zr 1 Mpc) 
that the form of the gamma function described above is stable for galaxies with r-band ab- 
solute magnitudes M a b s < —21. Our analysis of samples with even lower z max (the radius 
of the largest sphere is smaller than IOO/1" 1 Mpc) showed the stability of the presence and 
scale of the break at 20 — 25/i" 1 Mpc. 
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To assert that the distribution of galaxies is fractal up to 20h~ 1 Mpc, the power law must 
be observed in a considerably larger interval of scales (McCauley 2002), since inhomogeneous 
distributions of an arbitrary type can yield a power-law gamma function in a small interval 
of scales. At the same time, the homogeneity of the galaxy distribution on scales larger 
than lOO/i" 1 lie has been firmly established. The large-scale structure is characterized by 
the presence of inhomogeneities with an extent of more than lOO/i -1 Mpc. For example, the 
extent of the Sloan Great Wall structure is believed to be more than 300 Mpc. Therefore, 
the structuring of the Universe is not limited by the homogeneity scale, but, at the same 
time, the homogeneity of the distribution is consistent with the presence of structures on 
scales larger than the homogeneity scale (see, e.g., Gaite et al. 1999). 

Thus, the conclusions about the pattern of the distribution of galaxies, clusters, and 
superclusters reached by Tikhonov et al. (2000) are confirmed completely. 

5. THE LUMINOSITY DEPENDENCE OF CLUSTERING 

In this paper, we studied the dependence of the slope of the gamma function before 
the break on the luminosity range of the galaxies included in our volume-limited sample. 
To study the dependence of the correlation index 7 over a wide luminosity range, we chose 
the redshift boundaries z min = 0.03 and z max = 0.065. The upper limit on the absolute 
magnitude of the galaxies included in our volume-limited sample with these boundaries is 
Mu m = —19.58. We chose the Mu m range in which we analyzed the dependence, —19.6 

(14703 galaxies) 21.6 (1381 galaxies), in such a way that the sample was volume-limited 

and that a considerable number of galaxies remained in the sample (the number of objects 
in the sample decreases with decreasing Mn m ). In Fig. 2, 7 is plotted against Mi im . Since 
the samples drawn in this way overlap in bright objects, 7 varies smoothly. The dependence 
exhibits a break near Mu m = —20.2 followed by a more rapid increase in 7. 

To analyze the possible systematic trends related to the change in the number of objects 
and in the volume of space covered by the sample, we analyzed the dependence of 7 on the 
chosen far boundary (the redshift boundary z max )&t fixed Mu m = —21.5 (Fig.3) and the 
dependence of 7 on the longitude sample boundary \ max at \ m in = —47°, z max = 0.12, and 
Mn m = —21.5 (Fig.4). The index 7 increases as z max increases to 0.1. This may be explained 
by the fact that new significant structures fall into the sample. We can note a slight decrease 
in 7 with increasing \ max beginning from \ max = —10° as the volume and the number of 
galaxies increase. 

The dependence of the correlation index 7 on M abs in Fig. 5 was obtained using a sample 
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with the same geometrical boundaries as those in Fig.2, but with galaxies taken from nonover- 
lapping ranges of absolute magnitudes (the range is 0.2 wide). This increases the statistical 
significance of the features in the dependence. The number of objects smoothly decreases 
from 1861 in the range —19.6 < M^g < —19.8 to 555 at M abs < —22.0). The main feature of 
the dependence in Fig. 5 is that 7 increases significantly near the characteristic value of the 
luminosity function for SDSS galaxies, M* = —21.38 (at Hq = 65 km s _1 Mpc _1 )(Blanton et 
al. 2003), and above. 

6. CONCLUSIONS 

Thus, as the amount of data obtained as part of the SDSS project increases, it has 
become possible to directly test the cosmological principle. The homogeneity of the galaxy 
distribution has passed from a theoretical assumption confirmed by circumstantial evidence 
(isotropy of the cosmic microwave background, homogeneity of the sky distribution of radio 
and X-ray sources, etc.) to a strict observational fact, which automatically closes all the 
inhomogeneous models of cosmological evolution. The mean galaxy density is now a directly 
measurable quantity. 

Our results suggest that the slope before the break in the gamma function on a scale of 
~ 20/i _1 Mpc (the correlation index 7) depends in a complex way on the luminosity and the 
region of space in which the measurements are made. This casts doubt on the possibility of 
describing the distribution of galaxies by a single power law (with a single index 7) on short 
scales. 

The significant increase in the correlation index 7 with decreasing M abs after the break 
in the luminosity function for SDSS galaxies is indicative of a close relationship between the 
spatial distribution of galaxies and parameters of the luminosity function for these galaxies. 
The increase in 7 with luminosity must be reproduced in model calculations within the 
framework of the theory of galaxy formation. 
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Fig. 1. — Gamma function of a volume-limited sample of galaxies from the SDSS Main 
Galaxy Sample with z min = 0.05, z max = 0.16, M abs < —21.77, N = 18916 is the number 
of galaxies in the sample, R s = 132 Mpc, Ai is the break scale, and A 2 is the homogeneity- 
scale. 



-10- 




Fig. 2. — Correlation index 7 vs. limit on the absolute magnitude of the galaxies, M Hm , 
included in the sample with 0.03 < z < 0.065, -47° < A < 8°, 9° < rj < 36°. 
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Fig. 3. — Correlation index 7 as a function of the far sample boundary z max at z min = 0.03 
and M Hm = -21.5. 
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Fig. 4. — Correlation index 7 as a function of the longitude sample boundary X max at 
A miri = -47°, z max = 0.12 and M lim = -21.5. 




Fig. 5. — Correlation index 7 vs. M a b s range for nonoverlapping ranges of absolute magni- 
tudes at 0.03 < z < 0.065, -47° < A < 8°, 9° < 77 < 36°. The value of 7 corresponds to the 
upper boundary of the absolute magnitude interval and the next value corresponds to the 
lower boundary in the direction of decreasing M abs . 



